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/?-Cu2V207: a spin-^ honeycomb lattice system 

Alexander A. Tsirlinp Oleg JansonLO and Helge Rosneio 
Max Planck Institute for Chemical Physics of Solids, Nothmtzer Str. 40, 01187 Dresden, Germany 

We report on band structure calculations and a microscopic model of the low-dimensional mag- 
net /3-CU2V2O7. Magnetic properties of this compound can be described by a spin-i anisotropic 
honeycomb lattice model with the averaged coupling J\ = 60 — 66 K. The low symmetry of the 
crystal structure leads to two inequivalent couplings Ji and ,J[, but this weak spatial anisotropy 
does not affect the essential physics of the honeycomb spin lattice. The structural realization of the 
honeycomb lattice is highly non-trivial: the leading interactions J\ and J[ run via double bridges 
of VO4 tetrahedra between spatially separated Cu atoms, while the interactions between structural 
nearest neighbors are negligible. The non- negligible inter-plane coupling Jx — 15 K gives rise to the 
long-range magnetic ordering at Tjv — 26 K. Our model simulations improve the fit of the magnetic 
susceptibility data, compared to the previously assumed spin-chain models. Additionally, the sim- 
ulated ordering temperature of 27 K is in remarkable agreement with the experiment. Our study 
evaluates /?-Cu2V207 as the best available experimental realization of the spin-i Heisenberg model 
on the honeycomb lattice. We also provide an instructive comparison of different band structure 
codes and computational approaches to the evaluation of exchange couplings in magnetic insulators. 

PACS numbers: 75.30.Et, 75.10.Jm, 71.20.Ps, 75.50.Ee 
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I. INTRODUCTION 

Quantum magnetism is one of the most exciting and 
promising fields for exploring exotic ground states and 
unusual low-temperature properties. A variety of models 
and lattices lead to different regimes, ranging from sim- 
ple, collinear long-range order to intricate and essentially 
quantum ground statesii One of the relevant examples 
is the honeycomb (hexagonal) lattice that comprises a 
two-dimensional (2D) network of regular hexagons. The 
spin-^ Heisenberg model on the honeycomb lattice shows 
strong quantum fluctuations due to the low coordination 
number of 3^^— A weak interlayer coupling stabilizes the 
Neel ordering with the reduced sublattice magnetization 
of 0.54 fiB (compare to 1 fiB for a classical system) 4i"— 
Modifications of the model dramatically change its prop- 
erties. For example, frustrating next-nearest-neighbor 
couplings induce a spin-liquid or a valence-bond-solid 
ground statei^i"— The exactly solvable Kitaev model sets 
a specific arrangement of Ising-type interactions on the 
nearest-neighbor bonds of the honeycomb lattice and also 
leads to a spin-liquid ground stated that even sustains a 
partial disorder.^ Despite the diversity of intriguing the- 
oretical predictions, experimental studies remain scarce 
due to the lack of appropriate model materials. 

Recent studies proposed several high-spin honeycomb 
lattice compounds, based on Ni+^ (Ref. Il3 ) and Mn"'"'' 
(Ref.fl^. The iridates MalrOs (M = Li, Na) bear spin- 
2 and likely show anisotropic exchange interactions on 
the honeycomb lattice J^ In contrast, the isotropic inter- 
action regime is typically found in Cu'*'^ or V^** com- 
pounds, well to be accounted by a Heisenberg model. 
The apparently hexagonal arrangement of Cu atoms 
in Na3Cu2Sb06 is sometimes taken as an evidence of 
the honeycomb-lattice magnetism,— However, the struc- 
tural distortion and the orbital state of Cu make dif- 
ferent bonds of the lattice inequivalent and induce one- 



dimensional magnetic behavior.— A more appropriate 
system could be InCu2/3Vi/303 but it shows intrinsic 
disorder due to the intermixing of spin- 2 Cu'*'^ and non- 
magnetic V"*"^ cations. Although the Cu atoms are be- 
lieved to form small magnetic clusters on a honeycomb 
latticCj^iiiS, the influence of the non-magnetic part of this 
lattice is unclear and impedes the decisive comparison be- 
tween experimental results and theoretical predictions. 
Thus, structurally ordered and well-characterized spin- 
i honeycomb-lattice materials with essentially isotropic, 
Heisenberg exchange are still lacking. 



In low-dimensional magnets, the spatial arrangement 
of magnetic atoms is often deceptive, since the mag- 
netic couplings do not show a clear correlation with dis- 
tances between these atoms. In contrast, symmetries 
and overlaps of individual atomic orbitals play a cru- 
cial roleJ^"— In the following, we provide an instructive 
example that demonstrates the relevance of these char- 
acteristics for the microscopic magnetic model. We show 
that the /3-modification of CU2V2O7, previously consid- 
ered as a spin-chain or a spin-dimer compound, ^'^^^^ is 
the best available realization of the spin- 2 Heisenberg 
model on the honeycomb lattice. To arrive at this un- 
expected conclusion, we first examine the crystal struc- 
ture and analyze available experimental results fSec.llI|. 
Based on qualitative arguments, we rule out the previ- 
ously assumed spin-chain scenario and further evaluate 
the relevant exchange couplings along with the appro- 
priate spin model in a microscopic study (Sec. IIV[) . In 
Sec. |Vl we perform simulations for the proposed spin 
model and confirm our results by a direct comparison 
to the experimental data. We conclude our work with 
Sec. IVII that reviews the unusual implementation of the 
honeycomb spin lattice in /3-CU2V2O7 and provides an 
outlook for future experiments. 
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FIG. 1. (Color online) Crystal structure of /3-CU2V2O7 with copper and vanadium polyhedra colored in green (dark) and yellow 
(light), respectively. The left panel shows chains of edge-sharing CuOs pyramids running along [110] and [110]. The right panel 
presents a single layer of CUO4 plaquettes in the be plane, also indicated by a shaded area in the left panel. Individual exchange 
couplings are labeled by dark solid (Ji), dashed (Ji), light solid (Jd), and dotted (Ja) lines. The bottom part of the figure 
presents different representations of the Cu local environment: the CuOs pyramid (left) and the CUO4 distorted plaquette 
(right). The spin lattice is shown in Fig. [O] 



II. CRYSTAL STRUCTURE AND MAGNETIC 
PROPERTIES 



Experimental studies of /3-CU2V2O7 reveal low- 
dimensional magnetic behavior with predominantly an- 
tiferromagnetic (AFM) couplings (the Curie- Weiss tem- 
perature 6* ~ 82 K),^"* The temperature dependence of 
the magnetic susceptibility has a broad maximum at 
50 K, followed by the magnetic ordering transition at 
T/v = 26 K.^* The compound is controversially ascribed 
to different spin models. Two experimental studies ap- 
pHed the uniform-chain description ;^'^'^^ whereas Touai- 
her et al^ claimed a better susceptibility fit presuming 
the spin-dimer model. Further on, band structure cal- 
culations suggested two alternating couplings along the 
spin chains^S. Surprisingly, all the three proposals are 
quite different from the actual spin lattice. Below we 
demonstrate that empirical considerations in the experi- 
mental studies^i^i^^ do not take account of the peculiar 
electronic structure of /3-CU2V2O7, while the computa- 
tional study^^ failed to resolve all the relevant exchange 
interactions in the compound. 

The crystal structure of /3-CU2V2O7 is framework- like. 
Cu atoms have five-fold square-pyramidal coordination 
which is better described as "4-t-l" (see the bottom panel 
of Fig. [1]). Four oxygen atoms in the he plane reveal 
shorter bonds to Cu (1.93 - 1.95 A)^''' and form a CUO4 
plaquette, typical for Cu"*"^. The fifth oxygen atom shows 
a longer bond of 2.26 A which is roughly perpendicular 
to the plaquette along the a axis. The distance between 
the apical oxygen atom and the mean plane of the CUO4 
plaquette is about 2.2 A. Considering the Cu polyhedron 
as a pyramid, one finds edge-sharing connections to one 



neighbor along the h direction and to one neighbor along 
a, thus chains along [110] and [110] are formed (upper left 
panel of Fig. [T]). As the representation of the Cu environ- 
ment is reduced to the plaquette (right part of Fig. [1]) , 
the connection along a is broken, and we find CU2O6 
dimers of edge-sharing plaquettes (upper right panel of 
Fig. [T|) . Single VO4 tetrahedra link such dimers within 
the be planes, while the pyrovanadate [V2O7] groups join 
the resulting layers into a framework. 

Overall, the structure looks quite complex. A first 
guess on the possible magnetic interactions is to consider 
shortest Cu-Cu distances for the edge-sharing pyramids: 
2.95 A along h (the interaction Jd within the structural 
dimers) and 3.26 A along a (the interaction Ja). The 
resulting bonds form an alternating Ja — Jd chain pro- 
posed in Ref. [2^ (left panel of Fig. [1]). Other models 
assume the dimer limit {Jd 3> Ja)'^^ or the uniform chain 
limit Jd — Ja^^^ However, neither of these assumptions 
is correct due to the inappropriate choice of the leading 
couplings. The formation of a CUO4 plaquette is a clear 
signature of the magnetic (half-filled) orbital of d^2_y2 
symmetry, typical for Cu^^ with its 3(P electronic con- 
figuration. In /3-CU2V2O7, the plaquettes are slightly 
distorted. However, such a distortion has little effect on 
the crystal-field splitting, as confirmed by band structure 
calculations for /3-CU2V2O7 (Sec. lIVp and for other Cu+^ 
compounds. ^^ 

The magnetic orbital lies in the plane of the plaque- 
tte (be) and does not overlap with the orbitals of the 
fifth, apical oxygen atom. This feature, further confirmed 
by band structure calculations (Fig. 2]), allows us to re- 
duce the local environment of Cu to the CUO4 plaquette, 
despite the presence of five Cu-0 bonds in the crystal 



structure. The position of the ma gnet ic orbital makes 
the coupHng J^ negligible (see Ref. [2^ for a similar ex- 
ample). The coupling Jd is still allowed for the d^^^-y^ 
orbital. However, this coupling corresponds to the Cu- 
0-Cu angle of 98.7° implying sizable AFM and ferromag- 
netic (FM) contributions that might cancel each other.— 
These semi-empirical arguments are readily confirmed by 
our band structure calculations (Sec. IIVP and hint at rel- 
evant exchange couplings beyond the structural nearest 
neighbors. 

After resorting to the CUO4 plaquette description and 
considering the be plane as a potential magnetic layer, 
we find a remarkable similarity to the (VO)2P207 struc- 
ture (upper right panel of Fig. [IJ ^ This similarity is not 
surprising due to the similar (effective) A2B2X7 compo- 
sition. In fact, /3-CU2V2O7 and (VO)2P207 also share 
similar interpretation problems, because the symmetry of 
the magnetic orbital determines possible superexchange 
pathways and renders some of the short V-V contacts 
magnetically "inactive". Thus, the early misleading 
conjecture^ on the spin ladder physics in (VO)2P207 
was based on the incorrect assumption of the strong cou- 
pling via the V-O-V superexchange pathway, although 
this pathway lacks any magnetically active orbital (sim- 
ilar to Ja in /3-Cu2V207)J^ As one further employs the 
structural analogy between /3-CU2V2O7 and (VO)2P207, 
a strong AFM coupling J[ and alternating J( — Jd chains 
along the h direction could be expected (right panel of 
Fig. [T|). Yet, certain features of the electronic structure 
are different, thus even the J[ — Jd scenario of (VO)2P207 
has to be modified. Below, we show that the rearrange- 
ment of the cations (vanadium occupies the tetrahedral 
position and becomes non-magnetic, copper takes the 
vanadium position) strongly affects the electronic struc- 
ture. As a result, the spin system becomes 2D and attains 
the honeycomb-lattice geometry. 



III. METHODS 

To obtain a reliable microscopic model of /3-CU2V2O7, 
we perform extensive density functional theory (DFT) 
band structure calculations using the full-potential local- 
orbital scheme (FPLO9.00-33)ii2i We apply the local den- 
sity approximation (LDA) with the exchange-correlation 
potential by Perdew and Wang.i^ To cross-check the re- 
sults, we used three alternative computational schemes: 
i) the generalized gradient approximation (GGA)^ for 
the exchange-correlation potential within FPLO; ii) the 
Vienna ab initio simulation package (VASP)^ with 
the basis set of projected augmented waves;^ iii) the 
tight-binding scheme for linearized muffin-tin orbitals 
in atomic spheres approximation (TB-LMTO-ASA),'^^ 
where a different procedure for the evaluation of the ex- 
change couplings is implemented. The typical k mesh 
included 1098 points in the symmetry-irreducible part 
of the first Brillouin zone for the crystallographic unit 
cell (LDA calculation) and 288 points for the supercell 



(DFT-I-/7 calculations). We use the structural data3i 
from Ref. [33 as well as the data from Ref. JS^. The results 
for the two structures match very well, with the largest 
deviation below 10 % for one of the leading exchange 
couplings. 

The bare LDA approach does not provide a realistic 
description for the electronic ground state of transition 
metal compounds due to the underestimate of strong cor- 
relation effects in the 3d shell. Nevertheless, LDA results 
can be taken as a reliable input for the modeling. Follow- 
ing this idea, we extract the relevant LDA valence bands 
and fit them with a one-orbital tight-binding (TB) model 
using Wannier functions (WFs) centered on Cu sitesi^ 
The application of the WF technique leads to the unam- 
biguous evaluation of hopping parameters ti and provides 
a clear picture of the magnetic orbitals. Further on, we 
introduce the hopping parameters into a Hubbard model 
with the effective on-site Coulomb repulsion C/off. The 
Hubbard model is then reduced to the Heisenberg model 
to describe the low-lying excitations. This reduction is 
justified by the half-filling and by the strongly correlated 
{ti <^ C/off) regime. The parameters of the Heisenberg 
model are expressed as Atf/UcS and describe the AFM 
part jf^^^ of the exchange. The parameter UeS is fixed 
at 4.5 eV, according to previous studies i^'2^'^ 

The second approach to the evaluation of the ex- 
change couplings includes a mean-field treatment of cor- 
relation effects via the local spin density approximation 
(LSDA)-f [/ (or the related GGA+U) method. Two op- 
tions are possible: i) total energies for a set of coUinear 
spin configurations are mapped onto a classical Heisen- 
berg model (supercell approach); ii) the exchange in- 
tegrals are treated as second derivatives of the energy 
with respect to the rotation of the magnetic moments; 
such derivatives are calculated via matrix elements of the 
Green's function [Lichtenstein exchange integral proce- 
dure (LEIP)]4i The former approach is realized in FPLO 
and VASP, while the latter is implemented in the TB- 
LMTO-ASA code. The on-site exchange parameter of 
the LSDA+C/ method was fixed at Jsd = 1 eV. The 
choice of the Coulomb repulsion parameter U^d is fur- 
ther discussed in Sec. IIVBI 

The DFT-based microscopic model was challenged 
by the experimental magnetic susceptibility data and 
the magnetic ordering temperature, taken from Ref. |24| . 
The respective quantities as well as the high-field mag- 
netization curves, static structure factors, and spin 
correlations were computed via quantum Monte-Carlo 
(QMC) simulations with the loop algorithm^ or the di- 
rected loop algorithim in the stochastic series expansion 
representation;^ as implemented in the ALPS simulation 
package;-- The simulations are done for finite lattices 
with periodic boundary conditions. The typical lattice 
size was 16 x 16 (512 sites) for the 2D model (magnetic 
susceptibility, high-field magnetization, spin correlations, 
ordered moment) and up to 16 x 16 x 12 (6144 sites) 
for the three-dimensional (3D) model (magnetic ordering 
temperature) . The results obtained on lattices of differ- 
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FIG. 2. (Color online) Total and atomic-resolved LDA DOS 
for /3-CU2V2O7. The Fermi level is at zero energy. 



ent size ensured the lack of appreciable finite-size effects. 
To determine the ordered moment, we calculated static 
structure factors at a constant temperature T / Ji = 0.01 
and performed finite-size scaling, following Ref. |5|. 



IV. EVALUATION OF THE EXCHANGE 
COUPLINGS 

A. LDA-based model 

The main features of the LDA density of states (DOS, 
Fig. [2]) are typical for Cu"*"^ compounds j ^^i'^"'^^ Essen- 
tially, the valence bands are a mixture of Cu and O states 
with a small contribution of vanadium. These bands lie 
between —6.5 eV and 0.5 eV and account for the bond- 
ing between Cu and the five neighboring oxygen atoms 
in the CUO5 pyramid. The bands with the predominant 
vanadium contribution are found above 1.5 eV in agree- 
ment with the oxidation state of -1-5 for V atoms. The 
nearly isolated band complex at the Fermi level (Fig. [3]) 
is mostly formed by the magnetic Cu d^2_y2 orbital and 
the hybridizing 2po- orbitals of the four oxygen atoms 
comprising the CUO4 plaquette. The magnetic orbital 
lies in the plane of the plaquette and determines possi- 
ble superexchange pathways (see Sec. ^^. The energy 
spectrum is metallic due to the well-known shortcoming 
of LDA for strongly correlated systems. LSDA-|-[/ re- 
produces the insulating spectrum in agreement with the 
experimental identification of /3-CU2V2O7 as a magnetic 
insulator <22i2id^ 

The TB fit of the Cu dx2-y2 bands (Fig. [3]) yields three 
leading AFM interactions in the be plane: within the 
dimers (Jd), between the dimers along the b direction 
(J(), and between the dimers along [210] or [210] (Ji), 
see Table ID The nearest-neighbor interaction along a 
(Ja) is one of the leading magnetic couplings between the 
be planes, but it is much weaker than the next-nearest- 
neighbor in-plane couplings Ji and J(. This supports 
our empirical, qualitative considerations in Sec. |lTl Con- 



TABLE I. Cu-Cu distances (in A), leading hopping param- 
eters ti (in meV) of the TB model, and the resulting AFM 
contributions to the exchange integrals j/^^^ (in K). The ex- 
change integrals are derived as 4t|/[/eff with the effective on- 
site Coulomb repulsion C/cft = 4.5 eV. 
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FIG. 3. (Color online) LDA bands at the Fermi level (thin 
light lines) and the fit of the tight-binding model (thick 
dark lines). The notation of the k points is as follows: 
r(0,0,0), X(0.5,0,0), F(0, 0.5,0), Z(0, 0,0.5), r(0.5, 0,0.5), 
and A{0, 0.5, 0.5), and the coordinates are given along k^jky, 
and kz in units of the respective reciprocal lattice parameters. 



trary to the spin-chain scenario, we propose a 2D model 
with leading exchange couplings in the be plane. At first 
glance, the model approach would support the spin-dimer 
interpretation of Ref. [2y with the largest coupling J^. 
However, the nearly 90° superexchange regime of Jd im- 
plies a large FM contribution that strongly modifies the 
TB scenario (see Sec HVB]) .'^^ 

The couplings between the be planes are below 15 K. In 
addition to the nearest-neighbor interaction Ja, we find 
the long-range couplings J± with similar energy. Despite 
the short Cu-Cu distance, Ja is small, since the magnetic 
orbital of Cu does not overlap with the orbitals of the api- 
cal oxygen atom and makes the Cu-O-Cu superexchange 
impossible. The couplings J± run via the V2O7 groups 
along (approximately) [111] and [111] directions. Fur- 
ther hoppings in the TB model are below 25 meV, i.e., 
the respective J^^'^'s do not exceed 5 — 7 K and can be 
neglected within the minimum microscopic model. 

The calculated WFs (Fig. [J]) give a visual representa- 
tion of the orbitals contributing to the superexchange. 
Each WF includes the Cu d^^-y^ orbital and the p^ oxy- 
gen orbitals along with the smaller contribution of vana- 
dium orbitals. Since all the orbitals lie in the be plane. 




FIG. 4. (Color online) The Wannier function centered on the 
Cu site. 



the interplane couplings are relatively weak. Regarding 
the in-plane couplings, we note that oxygen contributions 
to WFs on the neighboring Cu sites show a nearly 90° 
arrangement for Jd- A similar arrangement is found for 
vanadium contributions with respect to Ji and J[. This 
causes a FM contribution to the exchange due to the 
Hund's coupling on the ligand site. 

Mazurenko et al.'^ developed a model for the ideal 90° 
superexchange geometry and calculated the FM part of 
the exchange: J^^ = —2/3'^JhNi, where /3 is the ligand 
contribution to the WF, Jh is Hund's coupling on the lig- 
and site, and Ni is the number of ligands where the WFs 
overlap (see also Ref . Wh . In our case, /Sq = 0.10 — 0.12, 
/3^ ~ 0.03, and Ni = 2. Assuming Jh = 1.5 eV^iil we 



find JJM 



850 K and Jf^^ ~ Jf ^ ~ -60 K. These 
values are rather overestimated (compare Tables HI and HIl 
jFM ^ _220 K, J™ ~ -10 K, and J'™ ~ -15 K), yet 
they identify a sizable FM contribution to Jd along with 
order of magnitude smaller FM contributions to Ji and 
J{. The overestimate is likely caused by the complex 
geometry of the system: the Cu-ligand distances are in- 
equivalent, whereas the angles exceed 90° (e.g., the Cu- 
0-Cu angle for Jd is 98.7°), thus significantly reducing 
the Hund's coupling. 



B. DFT+U 

To evaluate full exchange integrals via DFT+U calcu- 
lations, we consider the relevant in-plane couplings Jd, Ji, 
and J{ , as well as the two interplane couplings Ja and Jx . 
The results are collected in Fig. [5] and Table [ITl Most of 
the couplings are weakly dependent on the computational 
procedure, but the estimate of Jd is highly sensitive to 
the Ud value and to details of the DFT+U method. The 
investigation of the respective methodological problems 
constitutes the main part of this section. 

We start with the FPLO results shown in Fig. \5\ 
Different exchange-correlation potentials (LSDA-f C/ vs. 
GGA+U, upper panel) have little effect on the exchange 




FIG. 5. (Color online) Exchange couplings Ji = (2Ji +J[)/3 
(circles) and Jd (traingles) calculated with FPLO at dif- 
ferent values of Usd, the Coulomb repulsion parameter of 
DFT-l-f/. The upper panel compares LSDA-fC/ and GGA+U 
for the AMF double-counting correction scheme, while the 
bottom panel compares two double-counting corrections for 
the LSDA-I-C7 functional (see text). 



TABLE II. Exchange coupUngs calculated within LSDA-l-f/. 
The columns list individual exchange integrals (in K), the 
Coulomb repulsion parameter Uad of LSDA+f/ (in eV), the 
band structure code, and the double-counting correction. The 
solutions consistent with the experiment are marked with the 
bold typeface. 
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couplings. In contrast, the double-counting correction 
(DCC) scheme of DFT+U has a stronger effect on J's 
and especially influences Jd (bottom panel of the fig- 
ure). The DFT+U method adds an explicit mean- 
field energy term to account for the Coulomb repul- 
sion in the correlated shell. Then, a DCC is required. 
The DCC term subtracts a part of the repulsion en- 
ergy that is already contained in LSDA (GGA). Sev- 
eral DCC schemes are available, but their differences re- 
main weakly explored j "^^''*^ In the FPLO code, two DCC 
schemes are implemented: around-mean-field (AMF)^ 
and fully- localized- limit (FLL).'"'^ 

Band structure calculations for a family of copper ox- 
ides have rather firmly established the proper U^d value 
of 6 — 7 eV to reproduce the experimental exchange cou- 
plings within AMF.^^iiSii^iS Indeed, Usd = 6 eV in AMF 
leads to a reasonable scenario (first row of Table |TT|. The 



interactions Ji, J{, and J± are nearly unchanged com- 
pared to the TB estimates in Table HI This implies weak 
FM contributions in agreement with the long-range na- 
ture of these couplings. In contrast, the nearest-neighbor 
coupling Jd is reduced almost down to zero due to the 
superexchange pathway via the 98.7° bond angle. The 
other nearest-neighbor coupling Ja also has a FM contri- 
bution and becomes weakly FM. Since \Ja\ < 5 K, this 
interaction can even be omitted in a minimum micro- 
scopic model. As the U^d value is increased above 6 eV, 
the interactions Ji and J( are slightly reduced, while Jd 
becomes FM. 

Switching to FLL, we find a different scenario (Table ITTl 
and the bottom part of Fig. [5]). The couplings Ji and 
J[ are close to the AMF estimates. However, the short- 
range coupling Jd is about 200 K at U^d = 6 eV. As U^d is 
increased to 10 eV, Jd becomes relatively smaller, resem- 
bling the AMF result at U^d = 6 eV. Overall, we can reli- 
ably estimate Ji, J{ = 50 — 100 K, while Jd can be either: 
i) weak AFM (AMF, Usd = 6 eV or FLL, L/gd = 10 eV); 
ii) strong AFM (FLL, C/grf = 6 eV); iii) strong FM (AMF, 
Usd > 8 eV). In Sec. |Vl we carefully analyze the experi- 
mental data that unambiguously favor the first scenario. 
Unfortunately, this choice can not be made on purely 
computational (i.e., ab initio) grounds. The DFT+U 
calculations are typically compared to the experimen- 
tal band gap and the ordered magnetic moment (sub- 
lattice magnetization). The sublattice magnetization of 
/3-CU2V2O7 has not been reported and, anyway, should 
be subject to quantum effects (Sec. El) that lie beyond 
DFT+U and preclude the reliable comparison. The ex- 
perimental estimate of the band gap {Eg = 1.8 eV— ) 
is in good agreement to both AMF {Eg — 2.0 eV) and 
FLL {Eg = 1.75 eV) resuhs at Usd = 6 eV. As Um is 
increased to 10 eV, the band gap is somewhat overesti- 
mated {Eg — 2.3 eV in FLL). Thus, the computational 
arguments would rather prefer a "typical" Usd of 6 eV to 
describe the correlated electronic system of /?-Cu2V207. 
Yet, neither AMF nor FLL can be chosen unambiguously 
at this point. 

The DCC problem is not unique to the FFLO code, 
although in other codes it may be left out due to the lack 
of alternative DCC schemes implemented. In VASP, only 
the FLL and the related Dudarev's approach^ to the 
DCC are available. According to previous studiesj^i^ we 
select U^d = 7 eV and find a remarkable agreement with 
the FFLO FLL resuhs at Usd = 6 eV (see Table HU. The 
1 eV difference in Usd is a typical offset due to the differ- 
ent basis sets (the Coulomb repulsion potential is applied 
to the atomic 3d orbitals which are basis-dependent). 
Further on, we used the LMTO code and determined U^d 
from the constrained LDA procedure.^*' The resulting U^d 
of 9.5 cV is in agreement with previous reports^idliS and 
also leads to the typical FLL result with the large J^i^ 
Overall, different band structure codes converge to the 
same FLL solution for the exchange integrals. In this so- 
lution, Jd largely exceeds Ji and J{ and contradicts the 
experimental energy scale (see Sec.|Vl. To get a realistic 




FIG. 6. (Color online) Crystal structure of /3-CU2V2O7 and 
the sketch of the spin model with the leading couplings Ji 
(solid lines) and Ji (dashed lines). Dotted lines show Jd, the 
third-neighbor coupling on the honeycomb lattice. Open and 
shaded circles denote Cu atoms with opposite spins in the 
magnetically ordered ground state. 



solution with small Jd, one should cither switch to the 
AMF DCC or take a larger Usd value of about 10 eV for 
FFLO (and even larger values for VASP and LMTO). 

The problem of the DCC choice is not specific to /3- 
CU2V2O7, although in this compound it is probably most 
evident. Recently, we have shown that the AMF and FLL 
fiavors of LSDA-|-t/ lead to different results for several 
Cu"*"^ compounds i^iii^i^ Similar to /3-CU2V2O7, experi- 
mental data then favor the results at Uad = 6 — 7 eV for 
AMF and at t/sd = 9 - 10 eV for FLL within FPLO.-^ 
This empirical recipe helps to adjust the U^d value/the 
DCC scheme and to obtain realistic estimates of the 
exchange couplings. Still, a more general and compre- 
hensive theoretical study of this computational effect is 
highly desirable. 

The methodological part of the DFT+C/ study can be 
summarized as follows: 

i) The DCC scheme of DFT+U has influence on the 
exchange couplings and especially affects the short-range 
interactions; 

ii) The difference between the AMF and FLL flavors of 
DFT+[/ can be partly balanced by adjusting U^d- How- 
ever, the FLL results require unusually large U^d values 
that look unjustified with respect to the typical estimates 
(e.g., from constrained LDA or from the comparison to 
the experimental band gap). 

Switching back to the low-dimensional magnetism, we 
are able to construct a consistent microscopic picture, 
based on the TB analysis and the LSDA-I-C/ calcula- 
tions. We establish a 2D spin model for /3-CU2V2O7. 
The leading AFM interactions are Ji and J[ that act 
in the be plane. The non-frustrated interplane couplings 
J± amount to 10 — 15 K. The nearest-neighbor coupling 
Jd is somewhat ambiguous. However, the estimate of 
Jd — 200 K clearly contradicts the experimental energy 
scale. A thorough analysis in Sec. |V] shows that Jd is 
essentially a weak coupling, well below Ji and J{. The 
resulting Ji — J{ spin system reveals three bonds per site 
and represents the anisotropic honeycomb lattice (Fig.[B]). 
The coupling Jd between the structural nearest neighbors 



connects third neighbors on the spin lattice. Despite their 
complex structural arrangement, J± interactions can be 
considered as uniform interplane couplings. 



V. SIMULATIONS AND COMPARISON TO 
THE EXPERIMENT 

To challenge the proposed spin model by the exper- 
imental data and to investigate the properties of this 
model, we perform extensive QMC simulations. We start 
with the magnetic susceptibility*'" that was previously fit- 
ted with the expression for the uniform chain model.^ 
This fit yields J ~ 79 K and shows slight deviations from 
the experiment both at high temperatures and near the 
maximum (Fig. [7]). The model of the isotropic honey- 
comb lattice with Ji = 66 K produces a significantly 
better fit which is in remarkable agreement with the ex- 
periment down to T/v (see the difference plot in Fig. [T]). 
Below T/v, the susceptibility is strongly dependent on the 
direction of the applied field and can not be reproduced 
within the Heisenberg model that assumes isotropic ex- 
change. 

The low symmetry of the /3-CU2V2O7 structure leads 
to a weakly anisotropic honeycomb lattice with in- 
equivalent couplings Ji and J{. However, this spa- 
tial anisotropy does not affect the magnetic suscepti- 
bility. We were able to fit the data down to Tjv with 
Ji = (2Ji + J[)/3 = 60 - 66 K for J[/Ji = 0.5 - 1.3 
(the averaging includes two Ji and one J{ , because there 
are two Ji bonds and one J{ bond per site). This situa- 
tion resembles the anisotropic frustrated square lattice^ - 
where the moderate spatial anisotropy has little effect on 
the magnetic susceptibility. The fitted g value (1.9 — 2.2 
depending on the field direction, in reasonable agreement 
with the powder-averaged 17 ~ 2.1 from electron-spin- 
resonance measurements?'^) weakly depends on J[/Ji 
and can not be used to evaluate the degree of spatial 
anisotropy. Irrespective of the precise J'l/Ji ratio, the 
averaged coupling Ji is in excellent agreement with our 
computational estimates of 50 — 100 K for Ji and J[ (see 
Tables U and m. 

To resolve the ambiguity regarding the size of Jd, we 
extended the model and introduced the third-neighbor 
coupling on the honeycomb lattice (right panel of Fig. [S]). 
A small AFM Jd {Jd/ Ji — 0.2) does not lead to any vis- 
ible changes in the fit and slightly renormalizes Ji. The 
effect of the moderate interplane coupling J±^ — 0.2Ji 
is similar to that of Jd- However, an AFM intradimer 
coupling on the order of Ji {Jd = Ji) notably reduces 
the susceptibility maximum and leads to a poor fit of the 
data below 80 K (Fig.©. A larger Jd of 150 - 200 K (as 
suggested by the FLL calculations, see Table HB clearly 
contradicts the experimental energy scale. Finally, the 
sizable FM coupling Jd would strongly frustrate the spin 
lattice and impede the long-range ordering. Our simu- 
lations for the unfrustrated lattice show excellent agree- 
ment with the experimental ordering temperature (see 
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FIG. 7. (Color online) Upper panel: magnetic susceptibil- 
ity of /3-CU2V2O7 and the fits of the honeycomb-lattice (solid 
line), the uniform-chain (long-dashed line), and the honey- 
comb lattice with AFM Jd (short-dashed line) models. Ex- 
perimental data are taken from Ref. 24'. Difference curves are 
shown under the figure. Bottom panel: magnetization curves 
for the isotropic honeycomb lattice (solid line), anisotropic 
honeycomb lattice (dash-dotted line, J[/Ji = 0.5), and the 
uniform chain (dashed line) calculated at T / J\ — 0.03. 
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FIG. 8. (Color online) Binder ratio of the staggered magne- 
tization {mi)/{m^)'^ for the isotropic honeycomb lattice with 
the interlayer coupling J±/Ji = 0.2. Different symbols cor- 
respond to finite lattices of different size (note that the unit 
cell comprises two sites). The shaded bar shows the Neel 
temperature Tn/Ji — 0.41. 



below), hence the scenario of strong FM Jd is unhkely. 
We conclude that the coupling J^ is weak, while the lead- 
ing couplings are Ji and J[ forming an anisotropic hon- 
eycomb lattice. 

The magnetic ordering temperature T/y enables fur- 
ther comparison to the experimental data. To achieve 
the long-range ordering, we switch to a 3D model with 
an interplane coupling J^ ~ 13 K ~ 0.2Ji (cf. Tabic lll| . 
In a 3D model, the ordering transition is evidenced by 
a kink in the susceptibility curve or by a sharp increase 
in the staggered magnetization. To get a more accu- 
rate estimate, we calculate the Binder ratio {m'^)/{miy 
where rris is the staggered magnetization. The Binder 
ratio demonstrates a significant change upon the transi- 
tion (Fig. [5]). The intersection of the curves obtained for 
finite lattices of different size can be taken as a numer- 
ical estimate of Tn-^ We find Tn = 0.41 Ji ~ 27 K in 
remarkable agreement with the experimentally observed 
T/v = 26 Ki^ This comparison further supports our spin 
model. Similar to the magnetic susceptibility, the moder- 
ate spatial anisotropy {J[/Ji = 0.5 — 1.3) has little effect 
on Tn (at J(/Ji = 0.5, Tn/Ji = 0.37; at J{/Ji = 1.3, 
Tn/Ji = 0.40). 

The temperature dependence of the magnetic suscep- 
tibility favors the honeycomb-lattice description of /3- 
CU2V2O7 and yields the averaged coupling Ji. How- 
ever, the magnitude of the spatial anisotropy J[/Ji can 
not be precisely estimated. For the anisotropic frus- 
trated square lattice compounds, high- field magnetiza- 
tion measurements are the most suitable tool to evalu- 
ate the spatial anisotropy.^- Therefore, we also simulated 
magnetization curves of /3-CU2V2O7 (see the lower panel 
of Fig. [7]). The saturation field is slightly reduced for 
J'll Ji < 1 due to the decrease in Ji. The uniform-chain 
scenario would lead to an even lower saturation field of 
130 T, compared to 145 — 155 T for the honeycomb lat- 
tice. However, the easily accessible field region up to 
60 T shows negligible differences between the simulated 
curves. We conclude that high-field magnetization mea- 
surements could be a helpful tool to study the system, 
though only once the region around the saturation field 
(130- 155 T) is reached. 

Since the honeycomb spin lattice in /3-CU2V2O7 may 
be slightly anisotropic (0.5 < J(/Ji < 1.3), we study the 
infiuence of this spatial anisotropy on the ground state. 
The isotropic (J(/Ji = 1) lattice leads to the largest m 
of about 0.54 ^b- The spatial anisotropy reduces the di- 
mensionality of the system and, consequently, reduces m 
down to 0.47 fiB for J'l/Ji = 0.5 and down to 0.44 fj,B for 
J[/Ji = 1.5 (Fig. [HI). Although the magnitude of the ef- 
fect is 15 — 20 %, it may be hard to resolve experimentally 
due to the low absolute values of m and the high uncer- 
tainty of experimental estimates (typically, above 10 %). 
We also note that the interlayer coupling will rather in- 
crease the dimensionality leading to a slight increase in 
the m value. 

Spin correlations are a sensitive probe for the ground- 
state properties. Therefore, we calculate the expecta- 
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FIG. 9. (Color online) Finite-size scaling for the ordered mo- 
ment (sublattice magnetization, m) of the anisotropic honey- 
comb lattice model. L is the period of the finite lattice (cluster 
size) and Sq is the static structure factor at the propagation 
vector q. 



tion values {SqSb.) for R = 1 — 7 running along the Ji 
chains and in the perpendicular direction (see the in- 
set of Fig. [TUl) . For the isotropic (ideal) honeycomb lat- 
tice, both paths yield the same correlations (Fig. [TUl bold 
lines). The spatial anisotropy Ji/Ji ^ 1 slightly affects 
the correlations: for J{/Ji < 1, the enhancement of spin 
correlations along the Ji chains is accompanied by the 
weakening in the perpendicular direction. In contrast, 
the correlations between the chains become stronger for 
J'll Ji > 1. Although such a behavior is not surpris- 
ing and follows the simple physical picture of strong and 
weak bonds, the similarity of the curves corresponding 
to different J'l/Ji is remarkable. Fig. [10] evidences that 
the spin correlations strongly resemble the behavior of 
the isotropic honeycomb lattice, even for a sizable spatial 
anisotropy J(/Ji. This suggests that despite the spatial 
anisotropy, /3-CU2V2O7 is a surprisingly good realization 
of the spin- 1/2 Heisenberg model on the honeycomb lat- 
tice. 



VI. DISCUSSION AND CONCLUSIONS 

Based on a thorough microscopic study, we interpret 
the magnetic behavior of /3-CU2V2O7 within the spatially 
anisotropic honeycomb model. Recent theoretical reports 
show that the properties of this model depend on the 
magnitude of the spatial anisotropy.^"* For strong dimer 
anisotropy {J{/Ji > 1.85), a spin gap is opened, typical 
for a system of weakly coupled spin dimers. In contrast, 
a weak anisotropy leads to the honeycomb-lattice physics 
with the Neel AFM ordering. Our estimates of individual 
exchange couphngs suggest J'l/Ji = 0.5 — 1.3 and place 
/3-CU2V2O7 in the region of the Neel ordering, consistent 
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FIG. 10. (Color online) Spin correlations for different spa- 
tial anisotropics Ji/Ji compared to the isotropic honeycomb 
lattice (bold lines) along the J\ chains (top) and in the per- 
pendicular direction (bottom). The inset shows the respective 
paths; Ji and J( are depicted by solid and dashed lines, re- 
spectively. 



with the experiment i^ An AFM coupling J^ will further 
stabilize this ordering, while a FM Jd will induce mag- 
netic frustration. Our estimate of \Jd\/ Ji < 0.2 suggests 
that the effect of Jd should be relatively weak. Nev- 
ertheless, theoretical studies of the Ji — J{ — Jd model 
with the non-obvious ferromagnetic third-neighbor cou- 
pling Jd are desirable and could be relevant for related 
systems. 

The magnetic structure of /3-CU2V2O7 should bear fur- 
ther signatures of the underlying honeycomb spin lat- 
tice. The pronounced two-dimensionality and the low 



coordination number (three bonds per site) will lead to 
a strongly suppressed ordered moment of about 0.45 — 
0.55 fiB ,— well below the classical value. A similar ordered 
moment is found in green dioptase with its peculiar 3D 
spin lattice that comprises three bonds per site only.^- 
In /3-CU2V2O7, the spin arrangement should feature an- 
tiparallel ordering along Ji and J{ that further leads to 
an antiparallel ordering within the structural dimers (see 
Fig. ini). The expected propagation vector is k = (with 
respect to the atomic lattice), but the C-centering of the 
atomic structure has to be broken. In contrast, the AFM 
ordering within the Ja — Jd chains does not violate the 
C-centering. Thus, neutron diffraction experiments are 
a feasible test for the proposed 2D model of ^-Cu2V207. 
Inelastic neutron scattering should be able to resolve in- 
dividual couplings and to give a direct confirmation for 
the leading couplings Ji and J[. Additionally, the exci- 
tation spectrum of the spin-i honeycomb lattice can be 
studied. 

/3-CU2V2O7 gives an instructive example of the non- 
trivial magnetic interactions in transition-metal com- 
pounds. The symmetry of the magnetic orbital along 
with the Cu-O-Cu angle of 98.7° in the structural dimer 
disfavor sizable exchange couplings between the struc- 
tural nearest neighbors. Then, the long-range cou- 
plings come into play. The situation largely resembles 
VO(HPO4)-0.5H2O {Rei.M) or Cu2(P03)2CH2 (Ref.© 
where spin dimers do not coincide with the structural 
dimers. Another relevant example is BiCu2P06 with its 
rungs of the spin ladder running between the structural 
ribbonsj^i^ Such intricate implementation of the long- 
range superexchange couplings is an abundant well of 
surprises, as in the present study. /3-CU2V2O7 reveals 
a spin- 2 honeycomb lattice, despite lacking any appar- 
ent structural features of the honeycomb geometry. We 
note that the isostructural CU2P2O7 and CU2AS2O7 com- 
pounds could show similar features and deserve further 
investigation. 

In summary, we have shown that /3-CU2V2O7, previ- 
ously considered as a spin-chain compound, should be 
consistently described as a honeycomb-lattice system. 
The leading couplings Ji and J{ run via the non-magnetic 
VO4 tetrahedra, while the couplings between the struc- 
tural nearest neighbors are weak. The averaged cou- 
pling Ji amounts to 60 — 66 K. The spatial anisotropy 
is relatively small and has little influence on the ordered 
moment as well as on spin correlations. The interlayer 
coupling J± ~ 13 K leads to the Neel antiferromag- 
netic ordering at T^ = 25 — 30 K. We propose that 
/3-CU2V2O7 shows a relatively low sublattice magnetiza- 
tion of 0.45 — 0.55 /is, typical for the spin- 2 honeycomb 
lattice. Further tests of our model should include high- 
field magnetization measurements, neutron diffraction, 
and inelastic neutron scattering. We conclude that /3- 
CU2V2O7 is a noteworthy experimental realization of the 
spin-^ Heisenberg model on the honeycomb lattice. The 
convenient energy scale and the lack of disorder make 
it a promising system to challenge experiment and the- 
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ory and to improve our understanding of low-dimensional 
magnets. 
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